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E. Rubiola, M. Lenczner, P.-Y. Bourgeois, and F. Vemotte 


Abstract —This article Introduces the counter, a frequency 
counter — or a frequency-to-dlgltal converter, in a different 
jargon — based on the Linear Regression (LR) algorithm on 
time stamps. 

We discuss the noise of the electronics. We derive the statistical 
properties of the ft counter on rigorous mathematical basis, 
including the weighted measure and the freqnency response. We 
describe an implementation based on a SoC, under test in our 
laboratory, and we compare the f2 counter to the traditional 11 
and A counters. 

The LR exhibits optimum rejection of white phase noise, 
superior to that of the 11 and A counters. White noise is the 
major practical problem of wideband digital electronics, both in 
the Instrument internal circuits and in the fast processes which 
we may want to measure. 

The Ct counter finds a natural application in the measurement 
of the Parabolic Variance, described in the companion article 
arXiv: 1506.00687 [physics.data-an]. 


I. Introduction and State of the Art 

The frequency counter in an instrument that measure the 
frequency of the input signal versus a reference oscillator. 
Since frequency and time intervals are the most precisely 
measured physical quantity, and nowadays even fairly sophis¬ 
ticated counters fit in a small area of a chip, converting a 
physical quantity to a frequency is a preferred approach to 
the design of electronic instruments. The consequence is that 
the counter is now such a versatile and ubiquitous instrument 
that it creates applications rather than being developed for 
applications, just like the computer. 

The term ‘counter’ comes from the early instruments, where 
the Dekatron Q, a dedicated cold-cathode vacuum tube, 
was used to count the pulses of the input signal at in a 
reference time, say 0.1 s or 1 s. While the manufacturers 
of electrical instruments stick on the word ‘counter,’ the new 
terms ‘Time-to-Digital Converter’ (TDC) and ‘Frequency-to- 
Digital Converter’ are generally preferred in digital electronics 
(see for example 0, 0. 

The direct frequency counter has long time ago been re¬ 
placed by the classical reciprocal frequency counter, which 
measures the average period on a suitable interval by counting 
the pulses of the reference clock. The obvious advantage is 
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TABLE I 

Some commercial TDCs and Time Interval Analyzers. 


Device / Brand Resolution Type Reference 


Acam TDC-GPX 

10 ps 

Chip 

12 


Brilliant 

< 1 ps 

PCI/PXI Card 

13 


Guidetech 

1 ps 

PCI/PXI Card 

14 


K&K 

50 ps 

Special purpose 

15 



Keysight 53230A 

20 ps 

Lab instrument 

17 


Maxim MAX35101 

20 ps 

Chip 

18 


SPAD Lab TDC 

22 ps 

Chip 

19 


Texas THS788 

8 ps 

Chip 

20 



that the ‘1/n counts’ quantization uncertainty is limited by the 
clock frequency, instead of the arbitrary input frequency. Of 
course, the clock is set by design at the maximum frequency 
allowed for the technology. 

Higher resolution is obtained by measuring fractions of the 
clock period with a suitable interpolator. Simple and precise 
interpolators work only at fixed frequency. The most widely 
used techniques are described underneath. Surprisingly, all 
them are rather old and feature picosecond range resolution. 
The progress concerns the sampling rate, from kS/s or less 
in the early time to a few MS/s available now. See Q for a 
review, and Q for integrated electronics techniques. 

• The Nutt Interpolator 0. 0 makes use of the linear 
charge and discharge of a capacitor. 

• The Frequency Vernier is the electronic version of the 
‘Vernier caliper’ commonly used in the machine shop. A 
synchronized oscillator close to the clock frequency plays 
the role of the Vernier scale 0-00) 

• The Thermometer Code Interpolator uses a pipeline of 
small delay units and D-type flip-flops or latches. To our 
knowledge, it first appeared in the HP 5371A time interval 
analyzer implemented with discrete delay lines and sparse 
logic | fTT| . 

Table 1^ provide some examples of commercial products. 

In the classical reciprocal counter, the input signal is sam¬ 
pled only at the beginning and at the end of the measurement 
time T. There results a uniformly weighted averaging of 
frequency. In the presence of time jitter of variance tr^ (clas¬ 
sical variance), statistically independent a start and stop, the 
fluctuation of the measured fractional frequency has variance 
(Ty = 2cr^/T^. Increased resolution can be achieved with fast 
sampling and statistics. A sampling interval tq = t jm, to 1 
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interpolator 1 » 



t 

clock 1 II 1 rH 

~i V 

-^ 



time 


stamp^ 



II. Instrument Architecture and Noise 

The use of time stamps for sophisticated statistics is inspired 
to the ‘Picket Fence’ method introduced by C. Greenhall for 
the JPL time scale p0| , pT[ . Figure [T] shows a rather general 
block diagram. Following the signal paths, we expect that 
all the noise originating inside the instrument contains only 
white and flicker PM noise. The reason is that the instrument 
internal delay cannot diverge in the long run. Notice that 
the divergence of the flicker is only an academic issue p2) . 
Integrating over 20^0 decades of frequency, the rms delay 
exceeds by a mere 15-20 dB the flicker coefficient. Conversely 
the reference oscillator, which in a strict sense are not a part 
of the instrument, and the signal under test can include white 
FM and slower phenomena. The analysis of practical cases, 
underneath, shows that the instrument internal noise is chiefly 
white PM. 


Fig. 1. Time-stamp counter architecture. 


enables averaging on m highly overlapped and statistically 
independent measures. In this way, one can expect a variance 
(Ty oc There results a triangle weighted averaging 

of frequency. These two methods are referred to as 11 and 
A counters (or estimators) because of the graphical analogy 
of the Greek letter to the weight function f^, @ The n 
and A counters are related to the Allan variance AVAR j^ , 
| |24| and to the Modified Allan variance MVAR ||^-||Z7l. 
Frequency counters specialized for MVAR are available as a 
niche product, chiefly for research labs (for example |T5)). 

Having access to fast time stamping and to sufficient com¬ 
puting power on FPGA and SoC at low cost and acceptable 
complexity, we tackle the problem of the best estimator to be 
implemented in a frequency counter. 

The Linear Regression (LR) turns out to be the right answer 
because it provides the lowest-energy (or lowest-power) fit 
of a data set, which is the optimum approximation for white 
noise. The LR can be interpreted as a weight function applied 
to the measured frequency fluctuations. The shape of such 
weight function is parabolic (Fig. [^. We call the correspond¬ 
ing instrument ‘H counter,’ for the graphical analogy of the 
parabola with the Greek letter, and in the continuity of the H 
and A counters HD, HD. The fl estimator is similar to the 
A estimator, but exhibits higher rejection of the instrument 
noise, chiefly of white phase noise. This is important in the 
measurement of fast phenomena, where the noise bandwidth 
is necessarily high, so the white phase noise is integrated over 
the bandwidth. 

The idea of the LR for the estimation of frequency is not 
new ITT),®. However, these articles lack rigorous statistical 
analysis. Another modern use of the LR has been proposed 
independently by Benkler et al. |29| at the IFCS, where 
we gave our first presentation on the H counter and on our 
standpoint about the Parabolic Variance PVAR. 


A. Internal clock distribution 


The reference clock signal is distributed to the critical parts 
of the counters by appropriate circuits. As an example, we 
evaluate the the jitter of the Cyclone III FPGA by Altera p3) . 
A reason is that we have studied thoroughly the noise of this 
device p4| . Another reason is that it is representative of the 
class of mid sized FPGAs, and similar devices from other 
brands, chiefly Xilinx p5) , would give similar results. 

White phase noise is of the aliased (/?-type, described by 

= - 


where vq is the clock frequency, and k « 630 /rrad is an ex¬ 
perimental parameter of the component. Since the PM noise is 
sampled at the threshold crossings, i.e., at 21 ^ 0 . the bandwidth 
is equal to vq. Converting S^{f) into S^{f) = 4 ^S',^(/) 
and integrating on i^q, we get 


= 4^ 

thus Xrms = 2.8 ps at 100 MHz clock frequency. 

Flicker phase noise is of the x-type. We measured a value 
of 22 fs/vTlz at 1 Hz, and referred to 1 Hz bandwidth. 
Integrating over 12-15 decades, we And 115-130 fs rms, which 
is low as compared to white noise. 


B. Input trigger 

In most practical cases the noise is low enough to avoid 
multiple bounces at the threshold crossings, as described in 
I® In this condition, the rms time jitter is 

_ hj-ms 

Xrms = 

where 14ns is the rms fluctuation of the threshold, and SR is 
the slew rate (slope, not accounting for noise). 

The rms voltage results from white and flicker noise. The 
former is described as 14 = Cns/B, where e„ is the white 
noise in 1 Hz bandwidth, and B the noise bandwidth. The 
latter results from Vf = ae„ yJ\n{B/A), where ae„ is the 





























E. RUBIOLA &AL.THE Q. COUNTER, A FREQUENCY COUNTER BASED ON THE LINEAR REGRESSION JUNE 17, 2015 3 


flicker noise at 1 Hz and referred to 1 Hz bandwidth, and A 
is the low cutoff frequency. 

The golden rules for precision high-speed design suggest 

• White noise, e„ = 10 nV/V^Tz, including the input pro¬ 
tection circuits (without, it would be of 1-2 uV/VHz). 

• Flicker noise, a = 10... 31.6 (20-30 dB) at most, as a 
conservative estimate 

• the noise bandwidth is related to the maximum input 
(toggling) frequency i/max by i/^ax = 0.3 B. 

Let us consider a realistic example where = 1.2 GHz, 
e„ = 10 nV/-\/Hz, and a = 31.6 (30 dB). Accordingly, the 
input bandwidth is i? = 4 GHz. In turn the integrated white 
noise is Ku = 632 /iV rms. By contrast, flicker is Vf = 1.66 
/iV rms integrated over 12 decades (4 mHz to 4 GHz), and 
Vf = 1.86 p,V rms integrated over 15 decades (4 /iHz to 4 
GHz). White noise is clearly the dominant effect. 

C. Clock interpolator 

Commercially available counters exhibit single-shot fluctu¬ 
ation of 1-50 ps [Tab. [I] and References herein]. Since the 
interpolator is reset to its initial state after each use, it is 
sound to assume that the noise realizations are statistically 
independent, which is white noise. At a closer sight some 
memory between operations is possible, which shows up as 
some flicker noise. However, the same issues about bandwidth 
already discussed in this Section apply, and we expect that 
flicker is a minor noise effect as compared to white noise. 


D. Motivation for the Linear Regression 

The LR finds a straightforward application to the estimation 
of the frequency zz of a periodic phenomenon sin[^(f)] from its 
phase <p(t) using the trite relation u = ^ ^ At is well known 
that the LR provides the best estimate z> in the presence of 
white noise. The term ‘best’ means that the LR minimizes the 
squared residuals. This property matches the need of rejecting 
the white phase noise, which is the the main concern in 
counters, as we will see underneath. 


HI. The Linear Regression 
A real sinusoidal signal affected by noise can be written as 

v{t) = Vq sin[^(f)] 


where Vq is the amplitude, vq is the nominal frequency, 
and <p{t) is a phase that carries both the ‘ideal time’ and 
randomness. The randomness in <p{t) can be interpreted as 
either a phase fluctuation or a frequency fluctuation 


4>{t) = 


27rz/of + p{t) 

27rz/of + /(Azz)(f) dt 


PM noise representation 
FM noise representation 


Of course and (p{t) are allowed to exceed ±7r. 

We prefer to derive the properties of the LR using the 
normalized quantities x(t) and y(f), as in Fig. 


x(f) = t + x(f) (phase time) (1) 

y{t) = 1-1- y(f) (fractional frequency) (2) 



Fig. 2. Principle of the LR counter, and definition of often used variables. 


where 


x(f) = ip{t)l2TTVo (3) 

y(f) = x(f). (4) 


The quantity x(t) is the time carried by the real signal, which 
is equal to the ideal time t plus the random fluctuation x(f). 
The quantities x(f) and y(f) match the phase time fluctuation 
x{t) and the fractional frequency fluctuation y{t) used in the 
general literature |j^, p4| , with the choice of the font as the 
one and only difference. 

Most concepts are suitable to continuous and the discrete 
treatise with simplihed notation. For example, x = t + x maps 
into Xk = tk + Xfc for sampled data, and into x{t) = t + x{t) 
in the continuous case. The notation (.) and (.,.) stands for 
the average and the scalar product, defined as (x) = ^ 
and {x, y) = Xkyk for time series, where n is the number 
of terms in the sum, or as (x) = ^ f x{t)dt and {x,y) = 
/ x{t)y{t) dt, where T is the integration time in the continuous 
case. The span of the sum and the integral will be made precise 
in each case of application. The norm is dehned as ||a:|| = 
a/ {x, x). The mathematical expectation and the variance of 
random variables are denoted by E{ . } and V{ . }. 

The problem of the linear regression consists in identifying 
the optimum value y of the slope p (dummy variable used to 
avoid confusion with y) that minimizes the norm of the error 
X — r/t. The solution is the random variable 


y = 


(x-(x) ,t-(tj) 

Ilf-(f) IP 


(5) 


The choice of a reference system where the time sequence 
is centered at zero, i.e., (t) = 0, makes the treatise simpler, 
without loss of generality. Accordingly, the estimator y reads 


y = 



(6) 


These choice are equivalent because (t) = 0, so (x — (x), t) = 
(x,f) - (x) (l,f) = (x,t). 
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IV. Basic Statistical Properties 
A. Unbiased estimate 

The LR provides an unbiased estimate of the slope, i.e., 

E{y} = 1 (unbiased estimate), 

even without need that the noise samples (or values) are sta¬ 
tistically independent. This is seen by replacing the expression 
of the phase 


y = 


(x,f) (f + x,f) 


= 1 + 


(x,f) 


E{y} 


_1 I (E{x},f) 


V{y} = 


(estimator variance). 


This is seen by expanding the variance 

V{y} = E{(y-E{y})2} 

= E{f}-E{y}2. 


But 




+ 2{t, f)(x, t) + (x, tf} 




Since the random variables x^ are independent, then 

E{(x,t)“)=i: tkti E{xfcX<>} 
fc,r 

= '^tl E{x^} 




We conclude 


v{y} = i + ^<xx^-i = 


For even m = 2p, the proof starts with tk = kro for k € 
{-P, 

11*11' = E <1 = *-.'E = 

k——p k——p k—1 

2 2 

= (2p + 1) (p + 1) = (m + 2) (m + 1) 

^2 


mr 

~12 


for large m. 


B. Estimator variance 

Adding the hypothesis that samples x^ (or the values x(f)) 
are independent, then the estimator variance is 
^2 


A similar calculation can be carried out for odd m = 2p -1- 1 
with tk = ik+ i)ro for k € {-p - 1, ...,p}. 

D. Continuous input signal 

In the case of input signal continuous over a symmetric time 
interval (—§, §), we get 

12(x,<) 


y = l + 


r3 ■ 


Evidently {t) = 0, so the first equality results from ||f|p 
t^/12. Moreover, for a white noise x 

Indeed, replacing y by its expression. 

But E{(x, f)} = (E{x},f) = 0, so 

"^{y} = 

^r/2 |•T/2 

= 117114/ / ^MtHs)}tsdtds. 

I Til J-T/2J-r/2 


Since E{x(f)x(s)} = a'^S{t — s). 


^2 rT/2 


12 ■ 


and the conclusion follows. 

V. Weighted Average and Filtering 
We show that the Omega counter can be described as a 
weighted average or as a filter. The impulse response is derived 
from the weight function. 


C. Uniformly spaced time series 

In the case of a uniformly spaced data series with a constant 
step To and with r = mro, then for large m it holds that 


i + FM 


V{y) 


mr^ 

12a^ 


A. Weighted measure 

Let us consider an instrument that measures y(f) averaged 
on time interval of duration r. We express the estimate y(t) 
available at the time t as the weighted average 


y(t) = / y(s)'u;(s - t) ds, 
Jr 


( 8 ) 


rriT 


2 ■ 


(7) 


where w{t) is a suitable weight function which we will identify 
later in this Section, ruled by 

w{t) =0 for f ^ {—T, 0} (interval) 

/ w{t) dt = 1 (normalization). 

Jr 
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-t‘/ 2 rjl w{{) = + tI2) 



Fig. 3. Weight functions of the Q counter. 




Fig. 4. Impulse response functions of the counter. 


These rules expresses the fact that the integration time has 
duration t ending at the time t, and that the normalization 
yields a valid average. 

Eventually, (|^ can be seen as a scalar product, or as an 
operator on y, or as a weighted measure as regularly used in 
physics and engineering, or as a measure in the mathematical 
measure theory. 

We first identify a function w{t) which satisfies 

y(f) = / x{s)w{s — t) ds. (9) 

Js. 


Then, remembering that y(f) = x(f) and using the integration 
by part formula / dt = — J f'(t)g{t) dt, we find 


j{t) = — w(s) ds. 

J —OO 


( 10 ) 


Since we prefer the simplified expression 0* for the LR, 
we work in the (—§, §) interval using the centered functions 


Wc{t) = w{t + t/2) 


;(f) = w{t + t/2). 


Replacing w(t) with Wc{t), Eq. 0 gives the estimate at the 
time t + t/2 


rr/2 

y(r/2) = / x{t)wc(t)dt. 

J-t/2 

The LR (Eq. 0) applied over (—§,§) gives 


y(T/2) = 


(x,f) 12 r/2 


/-r/2 


x(f) t dt. 


Comparing GD and G^, we find 
12 

Wc{t) = ^t for f G (~ii §)> 0 elsewhere 


( 11 ) 


( 12 ) 


(13) 


and, using Eq. ( [T0) l, 
3 


Wc(t) = 


2t 


1 - 


At 


21 


forfG(—1,|), 0 elsewhere. 


These functions are plotted in Eig. 


B. Frequency response 

Eor the general experimentalist, the noise rejection prop¬ 
erties of the counter are best seen in the frequency domain. 
The fluctuations y(<) and x(f) are described in terms of their 
single-sided power spectral densities (PSD) Sy{f) and 
and the counter as a LTl (linear time-invariant) system which 
responds with its output fluctuation y(f). The counter output 
PSD is 


Sy{f) = \HU)?Sy{f) 

Syif) = \Hif)\^S,if), 

where |7T(/)p and |iT(/)p are the frequency response for 
frequency noise and phase noise, respectively. 

The LTI system theory teaches us that H{f) is the Eourier 
transform of the impulse response h{t) 

and so h{t) O H{f). 

Notice that a time shift f —>■ f ± | introduces a term 
in the Eourier transform, which vanishes in |JT(/)p and 
|iT(/)p. So, we can use the centered version, hc{t) O Tfc(/) 
and hc{t) -O- Hdf), of the above functions. 

We start from the time-domain response of the estimator, 
which results from the convolution integral 


y(t) = y(t) * h{t) 

= / y{s) h(t — s) ds. 

Jm 

A direct comparison of the above to 0 gives h{t) = w{—t), 
hence (Eig. 0 


h^t) = - 


Af 


for fG (—§,§), 0 elsewhere. 


and 

*- 62 

-^4/) =-(tt/t cos(7r/r) - sin(7r/r)) 

(7r/rcos(7r/r) - sin(7r/T))^ . 






























E. RUBIOLA &AL.THE Q. COUNTER, A FREQUENCY COUNTER BASED ON THE LINEAR REGRESSION JUNE 17, 2015 6 



f 


Fig. 5. Frequency response of the Q counter, plotted for r = 1 s.. |Ffc(/)P 
is the response to fractional-frequency noise Sy{f), while \Hc{f)\‘^ is the 
response to phase-time noise Sx{f)- 


In the same way, comparing the convolution integral 
y{t) = x(t) * h{t) 

= / x(5) h{t — s) ds 


to 


j 

(|^ yields h{t) = w{—t). There follows (Fig. 

12 

hc(t) = -for f e (—§, §), 0 elsewhere (14) 


and 

=-3^^ (tt/t cos(7r/r) - sin(7r/r)) 

(7r/rcos(7r/T) - sm(7r/r))^ . 

The frequency transfer functions |iT(/)p and |iT(/)p are 
plotted in Fig. 


VI. Hardware Techniques 

We describe a hardware implementation. Dropping the 
normalization, the LR estimates the frequency 


ET=oi(^k-{ek))itk-{t,)) 


ET=o (tk - 


(15) 


by fitting the phase data 6k = 4>k /27r expressed as the fraction 
of period. 

It is worth mentioning that timing indexes form a series 
of known natural terms such that the pre evaluation of {tk} 
and its average {tk) should be considered. The result of the 
average depends on m, the number of iterations. With m, some 
pre-calculated registers can be transferred into a RAM. 

Fig. 0 shows the block diagram. The algorithm takes two 
steps, represented as ‘stages.’ The reason is that both average 
and the samples over m are needed at the same time. As a 



Fig. 6. Hardware implementation. 


consequence, the first stage is mainly devoted to the evaluation 
of the average and the transfer of the data stream into an 
external dual-port RAM, the second stage to the calculation 
of the LR coefficients. 

First stage 

1) Data transfer: The incoming data flow is transferred into 
an external RAM at the interpolator sampling rate I/tq. In the 
presented design, I/tq is the effective sample rate of the data 
stream. 

2) Accumulation: Assuming {dk} where each 9k is coded 
in M bits, the data stream is firstly extended to 2M bits and 
shifted by a selected amount of gain g = log 2 (m). This way, 
when the accumulation process is done, there is no need to 
divide the result to obtain the average, by getting rid of the 
g lower bits. This process preserves the data integrity as the 
gain satisfies g ^ M. 

When the average is calculated and the first RAM full, an 
‘end’ of process signal and the result are propagated into the 
second stage. A new process cycle starts with the next data 
stream that will be stored in the second RAM. Such a cycle 
lasts obviously t = rriTo. The second stage is running at the 
same time, with the same duration r. Apart from an initial time 
lag T, this mechanism enable to calculate the LR coefficients 
at the rate of I/tq (so with no data loss): stage 2 is calculating 
on the old m data while the first calculates on the newer data. 

Second stage 

Here the averages are known and the samples collected, the 
calculation of the LR parameters are performed within this 
block. 

3) Process 1: Differences calculations. The samples are 
collected from the RAM and 9i — {6k) is calculated. In parallel, 
ti — (tk) are also computed, with (tk) provided by the CPU. 
The calculation of f — (tk) is cost effective compared to the 
RAM resource usage. All the data are latched. 

4) Process 2: frequency estimation. This process is an 
asynchronous accumulation J2n(^k — {dk)){tk — {tk))- The 
slope a is obtained by dividing the preceeding result by the 
constant denominator term provided by the CPU and stored 
into a register, and v is latched. Then, the estimated parameters 
are resized and propagated. 


































E. RUBIOLA &AL.THE n COUNTER, A FREQUENCY COUNTER BASED ON THE LINEAR REGRESSION JUNE 17, 2015 7 


TABLE II 


M 

m 

#RAM 

% 

#DSP48 

% 

32 

1024 

2/140 

1.4 

4/220 

1.8 

32 

215 

64/140 

45 

4/220 

1.8 

64 

215 

128/140 

91.4 

13/220 

5.9 


Oracle 

A testbench of the design requires different data sets to 
correctly predict its behavior. For we have implemented a C- 
based version code of the design in fixed point arithmetics, 
and equivalent to the GNU-Octave LinearRegression function 
or the Levenberg-Marquardt algorithm (as embedded into 
Gnuplot for example). 

A hardware implementation works at a rate of I/tq = 250 
MHz. 

The chip area occupation is only limited by the number 
of iterations m which induce high RAM usage with a linear 
dependence, m should be carefully chosen, depending on the 
required resolution/decimation rate and the design complexity. 
Excessive RAM usage may lead to long propagation time, thus 
to potential data corruption. 

For practical reasons, we also estimate the intercept point 
(b in Fig. with b = {(pk) — a (tfe)) leading to a more generic 
linear regression filter instantiation. Table [n| summarizes some 
implementation reports on the principal items that are BRAMs 
and DSP48, and their respective occupation for a Zynq-based 
system of the Zedboard type. 


VII. Discussion 


A. Noise Rejection 

We summarize the noise response of the counters in the 
presence of a data series {x^} uniformly spaced by tq, or 
equivalently sampled at the frequency = I/tq. The mea¬ 
surement time is T = mro. The noise samples associated 
to Xfc are statistically independent. 

The analysis of the H and A counters provided in this 
Section is based on | |2T| , | [22l , with a notation difference 
concerning the A counter. 

In pT) , p2) , the measurement time spans over 2t, and two 
contiguous measures are overlapped by t. This choice is driven 
by the application to the modified Allan variance, having the 
same response modtTy(T) = to a constant drift Dy. 

Oppositely, in this article use the same time interval t for 
all the counters, and we leave the two sample variances to the 
companion article 0. 

7) n counter: The estimated fractional frequency is 


^ ^7n—l Xq Xjji—i Xq 

y =-=-• 

mro r 

The associated variance is 

VS} = ^ ^ 

mTg 

independent of the sampling frequency. 


(16) 


2) A counter: In the measurement time r = mro we take 
m/2 almost-overlapped samples y^. measured by H counter a 
over a time t 12 


Yfc = 


^■mj2+k 
mro/2 


^mj2+k ^k 

^72 


The variance associated to each y^, is 


mrp /4 ’ 


The A estimator gives y = — X]fc=o 
tor variance is 

VS}=«^- 

m'^Tn 


' y^,. Thus, the estima- 


16cr; 




(17) 


3) n counter: The estimator variance 0 can be rewritten 
in terms of the measurement time r = mro and the sampling 
frequency i/g = I/tq as follows 


V{y} 


m^Tg 


12sT^ ' 


(18) 


4) Comparison: Comparing the estimator variance V{y} as 
given by Eq. ([TT]) and ( |T8| ) we notice that the H estimator 
features the poorest rejection to white PM noise, proportional 
to 1/r^. The white PM noise rejection follows the l/r^ law 
for both the A and H estimators, but the uj estimator is more 
favorable by a factor of 3/4, i.e., 1.25 dB. However small, 
this benefit comes at no expense in term of the complexity of 
precision electronics. 

The sampling frequency = 1 /tq can be equal to the input 
frequency izg or an integer fraction of, depending on the speed 
of the interpolator and on the data transfer rate. Values of 1^ 
MHz are found in commercial equipment. 


B. Data Decimation 

We address the question of decimating a stream {yfc(T)} of 
contiguous data averaged over t, in order to get a new stream 
{Vj (''”)} contiguous data averaged over nr, preserving the 
statistical properties of the estimator. 

7) H counter: Decimation is trivially done by averaging n 
contiguous data with zero dead time 

- n—1 

yjnr) = - ^yfe(T). 
k=0 

2) A counter: Decimation can be done only in power of 2 
(n = 2^), provided that contiguous samples are overlapped 
by exactly r/2. In this case, we apply recursively the formula 

y(2r) = ^y-i(T) + ^yo(^) + ^yiM- 

Notice that the measurement time 2r is given by the overlap 
between samples. 

3) counter: An exact decimation formula, as for the 
H and A counters, does not exist. The reason is that the 
exact parabolic shape has constant second derivative. This is 
incompatible with the presence of the edges at ±r/2 when 
more than one shifted instances of the w{t) function are 
overlapped. 
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C. Application to the two-sample variances 

The general form of the two-sample variance reads 

V{y} = iE{(yi-y2)"} (19) 

Using a stream of ‘contiguous’ data {y^} measured with a If, 
a A or a n estimator, Eq. ( [T9| , gives the AVAR, the MVAR or 
the PVAR. Notice that in the case of the MVAR, ‘contiguous’ 
means 50% overlapped. 

VIII. Conclusions 

The same architecture (Fig. is suitable to the If, A and U 
estimators, just by replacing the algorithm that processes the 
phase data. A smart design can implement the three estimators 
with a some digital-hardware overhead. 

The If estimator is the poorest, and mentioned only for 
completeness. Nonetheless, it is still on the stage when ex¬ 
treme simplicity is required, or in niche applications like the 
measurement of the Allan variance. 

In the presence of the same amount of white noise, the U 
estimator is similar to the A estimator, but features smaller 
variance. The difference is a factor of 3/4 (1.25 dB). 

The A estimator is a good choice when a clean decimation 
of the output data is mandatory. By contrast, the U estimator 
is superior for the detection of noise phenomena | |J7| 

However may object that a noise reduction of 1.25 dB is 
not a big deal, the value our analysis is in showing that there 
is no room for further improvement. 

If complexity is the major issue, the A estimator fits, at the 
cost of 1.25 dB higher noise. Oppositely, if the ultimate noise 
performance in the presence of wideband noise is an issue, 
the H estimator is the right answer. The complexity of the LR 
is likely to be affordable on FPGA or SoC electronics. 
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